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We investigate controlled collisions between trapped but separated ultracold atoms. The inter- 
action between atoms is treated self-consistently using an energy-dependent delta-function pseu- 
dopotential model, whose validity we establish. At a critical separation, a "trap-induced shape 
resonance" between a molecular bound states and a vibrational eigenstate of the trap can occur. 
This resonance leads to an avoided crossing in the eigenspectrum as a function of separation. We 
investigate how this new resonance can be employed for quantum control. 

PACS numbers: 34.50.-s 03.67.Lx 32.80.Pj 34.90.+q 



The ability to arbitrarily manipulate the quantum 
state of a many-body ensemble represents the ultimate 
control of a physical system. This task has steadily ad- 
vanced in atomic-molecular-optical systems with tremen- 
dous progress in cooling and trapping technology. This 
has led to the creation of Bose-Einstein condensates 
(BEC) and Fermi degenerate gases, and the explorations 
of new forms of matter and mesoscopic quantum states 

Heviously accessible only in condensed matter systems 
. The addition of engineered traps, such as optical lat- 
tices 01 and other optical Q and magnetic jj] microtraps 
provides a new knob with which to control the quantum 
state. A dramatic example of many-body control in lat- 
tices was demonstrated through the observation of a su- 
perfluid to Mott insulator quantum phase transition and 
the collapse and revival of the mean field coherence Q . 

The standard approach to modelling and designing co- 
herent states of matter, such as occur in a quantum phase 
transition, has its foundations in condensed matter the- 
ory, where one considers solutions to the entire many- 
body Hamiltonian. An alternative viewpoint arises from 
a fundamental theorem of quantum information theory 
0: an arbitrary state of a many-body system can be 
reached entirely through operations on single bodies and 
pairwise interactions. This provides a direct approach 
to engineering mesoscopic states through the application 
of a "quantum circuit" 0- Moreover, one requires only 
a single two-body interaction (e.g. CPHASE or CNOT 
gate) that entangles the "particles" to contribute to a 
universal set of quantum logic gates. 

In the context of ultracold neutral atoms, whereas ma- 
nipulating the quantum state of an individual atom is 
a very mature technique, arbitrary unitary mapping of 
a two-atom system has not yet been achieved. Neu- 
trals, by their very nature, do not strongly couple to 
anything. This may be an advantage for avoiding noise, 
but it implies that the two-body interaction will gen- 
erally require close overlap of the atomic wavepackets. 



By bringing two atoms within the same well of a tightly 
confining microtrap, one can achieve this strong coupling 
while remaining in the electronic ground state. Propos- 
als for two-atom control in such a geometry have been 
considered using ground state s-wave collisions Fes- 
hbach resonances and laser induced Raman transi- 
tions At such close range, the atoms lose their 
individual identities and instead must be described as 
a molecular dimer which generally does not respect the 
atomic symmetries. This constrains the possible encod- 
ings of quantum information such that two-body logic 
gates can be performed within a well-defined "logical ba- 
sis" . This constraint can be overcome by placing the 
particles in distinguishable locations where the atomic 
quantum numbers are conserved asymptotically. Under 
typical conditions such separated atoms would generally 
encounter very weak interactions. The coupling between 
atoms can be dramatically increased, however, when a 
resonance of the two-body system is excited, resulting 
in long-range interactions. An example of this is induced 
electric dipole-di pole int eractions associated with excited 
electronic states 
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In this letter we describe the physics at the founda- 
tion of these protocols by considering ultracold collisions 
between trapped but separated atoms. Our studies show 
that the trapping potential can lead to new resonances 
not found in free space. Unlike standard atomic colli- 
sions, here the relative motion of the atoms is quantized 
by the trap, making its description intrinsic to the pro- 
cess. Resonances can then occur between eigenstates of 
the trap and molecular bound states, allowing us to over- 
come the generally very weak interactions associated with 
the van der Waals potential. These "trap induced shape 
resonances" (TISR) can be substantial but are not ac- 
counted for within perturbation theory, that was applied 
in previous proposals for quantum logic via cold collisions 
The TISR provides a new tool for molecular dimer 
control (e.g., the production of cold molecules) and the 
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FIG. 1: Sum of the harmonic trapping potential and chemical 
binding potential (gray line) in the relative coordinate r for 
zero trap separation (a) and larger trap separation Az ^ 
(b). The molecular bound state at Eb and trap eigenstate at 
-Etrap can become resonant at a critical separation Azrcs. 



design of two-atom quantum logic gates. 

Our model system consists of two trapped and sepa- 
rated atoms that interact through the molecular poten- 
tial Mnt, described by a Hamiltonian, 



2m 

+Mnt(ri - r2) 
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where Az is the separation of the traps (chosen in the z- 
direction). The trapping potential Vt for the two atoms 
could be, for example, the state-dependent trap of a 
three-dimensional optical lattice potential. In this system 
Az can be continously controlled by the angle between 
the polarization vectors of the counter propagating laser 
beams 0, Q| • We assume atoms are well- localized near a 
potential minimum that is approximated as isotropic and 
harmonic with frequency uj. The two-atom Hamiltonian 
then separates into one for the center-of-mass moving in 
an isotropic harmonic potential, and one for relative co- 
ordinate dynamics governed by 



^ + iA^c.2|r-Az|VM„t(r) 



(2) 



The reduced mass fj, = m/2 moves under the combined 
effects of a harmonic trap centered at Az and a central 
interatomic potential (Fig.^. 

1 /2 

The harmonic trap is characterized by zq = {h/ fioj) , 
the width of the trap ground wave function, while the 
interatomic potential has a much shorter range a:o; for 
ground state van der Waals interaction Vint{r) — —C^/r^ 
at large interatomic separation, xo — (2/iC6/fi.^)^^^/2 
[l^ . Since for alkali atoms in typical optical lattices 
xq zq, distortion of the interatomic potential due to 
the harmonic trap can be neglected. Furthermore, for 
the ultra-cold collisions under consideration here, s-wave 
scattering dominates the collision. Under these condi- 
tions, the interatomic interaction Vint(r) can be replaced 
by a zero-range effective-scattering length pseudopoten- 
tial Illllllia 



ji or 



Here Ek = (/ifc)^/(2/z) is the kinetic energy of relative 
motion for two atoms in an asymptotic scattering state 
with momentum hk. It determines an "energy-dependent 
scattering length," 
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where So{Ek) is the s-wave collisional phase shift. How- 
ever, unlike in a traditional collision where the two atoms 
are asymptotically free, in this system the atoms are 
trapped. The eigenvalues of the system must thus be 
solved self- consistently 0|. To this end, the eigenspec- 
trum of the system is first calculated as a function of 
the (energy-independent) scattering length. Second, the 
effective-scattering length is calculated as a function of 
kinetic energy Ek for two untrapped atoms. The self- 
consistent energy eigenvalues are then found numerically 
as the simultaneous solutions of these functions. This 
method has been shown to yield the correct scattering 
behavior and energy spectrum for two atoms in a tight 
harmonic trap |0| for zq 3> Xq. 

Full characterization of the interatomic potential by a 
pseudopotential requires not only accurately reproduc- 
ing the scattering behavior but also the molecular bound 
state spectrum. This is particularly important for the 
TISR described below. We accomplish this by analytic 
continuation of the effective scattering length Eq. to 
negative kinetic energies Ek = —h?K'^/{2jjL) according 
to UesiEK) = —tAnh.%So{EK) / with real and positive 
K. The self-consistent solutions for these negative ener- 
gies then accurately reproduce the entire s-wave bound 
state spectrum as can be easily understood as follows. 
Suppose there is a bound state of the actual potential 
Vint at Eb = —{hKb)'^/{2ii). The corresponding pole of 
the S'-matrix implies that tanhi^o(£'6) 1, and hence 
o-csiEb) = l/Kfc. Now, a pseudopotential with scatter- 
ing length a > possesses exactly one bound state at 
Es = —{h?/{2^a'^). By setting a = aes{Eb) we recover 
the exact bound state, Es = Eb- Note that without a 
self-consistent solution one cannot choose a pseudopoten- 
tial which simultaneously matches the scattering length 
and last bound state of the true potential. Moreover, al- 
though the pseudopotential can have no more that one 
bound state, this self-consistent solution via the energy- 
dependent scattering length captures exactly all of the s- 
wave bound states, since they are obtained from Ooff (-E/f ) 
and therefore the S'-matrix. 

With this model we calculate the energy spectrum of 
Eq. (O for two trapped but separated atoms. First con- 
sider the spectrum for a given scattering length. We 
represent the Hamiltonian for arbitrary Az in the ba- 
sis corresponding to the solutions with Az = and a 
fixed scattering length a (i.e. not the self-consistent so- 
lution). This basis, derived by Busch et al. 0, consists 
of 3D-harmonic-oscillator-like solutions. The pseudopo- 
tential couples only s-waves, resulting in irregular I — 
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FIG. 2: Energy spectra as a function of separation Az be- 
tween traps for a negative (a) or positive (b) fixed scattering 
length a. The results of perturbation theory are shown as 
dashed lines. For a > one sees the parabolic energy shift 
of the molecular bound state due to the harmonic trapping 
potential and the avoided crossings associated with the TISR. 

radial waves with singularities at the origin. These so- 
lutions also include the pseudopotential bound state at 
negative energy. The I > 1 wave functions are the regular 
3D-harmonic oscillator wave functions. The trap poten- 
tial, proportional to |r Az|^ = — Azrcos0 + Az^ 
is axially symmetric and dipolar, thereby preserving the 
magnetic quantum number of relative motion and cou- 
pling the partial waves I to I ± 1. The Hamiltonian is 
diagonalized in this basis at each trap separation Az. 
The resulting energy spectrum for two atoms is shown in 
Fig. 121 for both positive and negative scattering lengths, 
a = +0.5zo and a = — 0.5zo, as a function of Az. The 
results of first-order perturbation theory are also shown 
for comparison. For Az ^ zq, we recover the expected 
unperturbed 3D harmonic oscillator eigenenergies. As 
the separation between traps is decreased, perturbation 
theory predicts a negative (a < 0) or positive (a > 0) en- 
ergy shift to the ground state. The expected behavior is 
seen for a = — O.Szq, but an unexpected solution is seen 
for a = +0.5zq. 

The results for the positive scattering length are ex- 
plained as follows. For large positive a there is a molecu- 
lar bound state close to dissociation. As Az is increased, 
the interatomic potential, located at very small internu- 
clear distances, is pushed up in energy by fj.to'^Az^ /2 due 
to the parabolic trapping potential in Eq. Q. That is, 
in order for the separated atoms to collide, they must 
overcome the potential barrier created by the trap (see 
Fig. nj . When the molecular bound state becomes reso- 
nant with the lowest trap eigenstate, an avoided crossing 
occurs in the energy spectrum(see Fig. |21l. As the sep- 
aration is increased even further, the molecular bound 
state becomes resonant with higher-lying trap states and 
more avoided crossings occur. This is a new "shape res- 
onance" for s-wave collisions in which the trap barrier 
plays the role of the centrifugal barrier in a standard free- 
space shape resonance for higher partial waves. Analo- 
gous, Feshbach-like, "confinement induced resonances" 
have recently been predicted for ID and 2D trapped 



bose eases and for delocalized states in 3D-optical lat- 
tices 

The separation at which the lowest resonance occurs, 
Azrcs, is easily estimated by equating the sum of the 
molecular binding energy and trapping potential at the 
origin, Ei, + /xw^Az^/2, to the vibrational ground state 
energy of the oscillator, 3huj/2, yielding Azies/^o = 
^/S + Zq /a^ . The location and gap of the avoided cross- 
ing depends strongly on the molecular binding energy. 
For a deeply bound state, but still positive scattering 
length, corresponding to < a ^ zq, the resonance oc- 
curs at much larger separations and with an exponen- 
tially small energy gap. This corresponds to the small 
probability for the atoms to tunnel from the trap into 
the chemical bin ding potential. Using a standard vari- 
ational approach |20| based on symmetric and antisym- 
metric combinations of the Az = bound state |0| and 
the trap ground state, we find that for < a < 0.2zo 
the gap is smaller than 10~^huj. For a zq the energy 
gap asymptotes to a large value, Ai?niax = 0.5640?iaj. 
The shape resonance can therefore be easily observed for 
large positive scattering lengths, where the bound state 
would be close to dissociation. This can be achieved in 
tight traps, where the scattering length is on the order 
of zo and the energy gap approaches a significant frac- 
tion of huj (see Fig.O. For example, in an optical lattice 
of ^^^Cs atoms the large scattering lengths of 280ao to 
2400ao are comparable to typical trap sizes in an optical 
lattice, corresponding to a modest Lamb-Dicke parame- 
ter 7] = k^zQ = 0.1. A substantial TISR will result. 

To obtain a more accurate spectrum in the case of 
trapped alkali atoms, we must account for the energy- 
dependence of the scattering length in the self-consistent 
model described above. As a test case we consider 
the simplest possible interatomic potential - a step- 
potential of radius R and depth Vq with a single s-wave 
bound state. The s-wave phase shift is given explic- 
itly by So (EktVo) = arctan (fc tan(gi?) /g) — kR, where 
q = ^/2fj.(Ei^ + Vn)/h?^I\. The energy-dependent scat- 
tering length is evaluated using Eq.l^J as a function of 
Ek- The relative kinetic energy of the colliding atoms 
is given by the total energy eigenvalue E minus the trap 
potential at the origin. 

Figure 13 shows the self-consistent energy spectrum as 
a function of well separation Az. These approximate 
eigenvalues are compared with the exact solution for the 
step-potential (Vq = 36.79huj and R = 0.2zo) plus har- 
monic potential, calculated numerically. We accomplish 
this by expanding the total Hamiltonian in isotropic 3D- 
harmonic oscillator wavefunctions and diagonalizing the 
matrix. Figure |2| also shows a plot of the constant scat- 
tering length approximation, using the zero-energy scat- 
tering length a = acff(O). As expected, this approxi- 
mation fails to capture the correct bound state energy 
and therefore the correct location of the shape reso- 
nance. In contrast, the self-consistent solution using the 
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FIG. 3: Comparison between the energy spectrum of the test 
step-potential and that of the pseudopotential approximation. 
The two lowest energy curves are shown for the step-potential 
(solid), the pseudopotential with an energy dependent scat- 
tering length flcff (circles) and constant a (dashed). The inset 
shows the energy gap Ai5 for these two lowest levels at res- 
onance as a function of scattering length a calculated from 
energy spectra at the different scattering lengths. A varia- 
tional estimate of the energy gap is shown as the dashed line. 

energy-dependent pseudopotential shows excellent agree- 
ment with the exact calculation, even for a well that was 
chosen to have a fairly long range, R = 0.2 2:0. The agree- 
ment only breaks down when the range of the potential 
becomes on the order of trap size, R > 0.5zo. 

Having confirmed the validity of the self-consistent 
pseudopotential model, we can use the calculated en- 
ergy spectrum to design two-atom quantum logic gates. 
The use of ultracold collisions was initially considered by 
Jaksch et aZ.Q in the non-resonant case, using pertur- 
bation theory. Our analysis shows that in principle, for 
positive scattering lengths, resonances will occur at some 
atomic separation, and perturbation theory will break 
down. The resulting avoided crossing in the energy spec- 
trum must be properly accounted for. This is particularly 
true for atoms with very large scattering lengths, such as 
^^■^Cs. The TISR opens the door to new protocols for 
entangling two-atom logic gates with separated atoms. 
For example, a 27r Rabi oscillation between the trapped 
atoms and an auxiliary molecular bound state leads to 
a phase shift of —1 on the two-atom state. If the ac- 
quired phase shift occurs only for one logical encoding of 
the atoms, the resulting unitary transformation is the so 
called "CPHASE" two-qubit logic gate 0. An optimal 
regime for operation of this protocol is where zq ^ a. In 
this regime the energy gap (FigOinset) asymptotes to its 
maximum value, minimizing the dependence of the coop- 
erative phase shift on the precise value of the trap width 
00, and hence reducing errors due to trap-laser intensity 
fluctuations. 

More generally, beyond quantum logic, the TISR pro- 
vides a new avenue for spectroscopy and coherent con- 



trol of ultracold molecular dimers. Like magnetic Fesh- 
bach resonances, these shape resonances can provide new 
ultra-high precision spectroscopic data on the molecular 
potential |2l| , and the production of cold molecules tun- 
able by the trap parameters. A complete characterization 
of these protocols requires a generalization of our model, 
including the full spin-dependent nature of the collision 
process via the hyperfine and exchange interactions. We 
plan to address this in future work using multichannel 
Born-Oppenheimer potentials approximated by energy 
dependent pseudopotentials. 
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